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The nonperturbative QCD vacuum at finite temperature and a finite baryon den- 
sity in an external magnetic field is studied. Equations relating nonperturbative 
condensates to the thermodynamic pressure for T ^ 0, fig ^ and H ^ are ob- 



o. 

psj I tained, and low-energy theorems are derived. A bulk viscosity ^(T, ^u, H) is expressed 



in terms of basic thermodynamical quantities describing the quark-gluon matter at 
r 7^ 0, ^ 0, and H ^ 0. Various limiting cases are also considered. 
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1. INTRODUCTION 



Quantum chromodynamics (QCD) as a theory that describes strong interaction physics is 
still vigorously developing. Over the past few decades, particular attention has been drawn 
to the investigation of the behavior of strongly interacting matter under the effect of various 



^ i external fields. In the real world, these are primarily temperature and the baryon density. 

Interest in the behavior of matter under extreme conditions (high temperature is the one 

starting from a characteristic QCD scale, T ~ 200 MeV, and high baryon densities starts 

^ I from n > tiq ^ 0.17 fm~^, where no - is a normal nuclear density) is motivated primarily 

H ! 

by the experiments studying heavy-ion collisions. Due to this, one can expect that such 
experiments probe densities and temperature at which a phase transition to quark-gluon 
plasma, which is a new state of strongly interacting matter, is possible. 

In recent years, the phase structure of the vacuum in an external magnetic field H has 
become yet another important object of investigations. It was shown recently that magnetic 
fields of strength in the range eH ~ 10 -f- lO'' MeV^ can be generated in heavy-ion collisions 
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[1-3]. Such field can lead to observable phenomena (so-called "chiral magnetic effect") [1, 4- 
7] in experiments at the Relativistic Heavy Ion Collider (RHIC) and at the Large Hadron 
Collider (LHC). Magnetic fields of the order of eH ~ ^%cd or even stronger could exist in 
the early Universe. Such magnetic field strengths can lead to new interesting phenomena 
accompanying the QCD phase transition [8-32]. 

In [33], the magnetic- field dependence of the quark condensate was studied on the basis 
of the Nambu-Jona-Lasinio model. In QCD, the one-loop result for the H dependence of 
{qq) was obtained in [34]. In both cases, the condensate was found to grow with increasing 
H. It follows that a naive analogy with superconductivity theory, where a magnetic field 
breaks the condensate of Cooper pairs, is inapplicable here. Strictly speaking, the behavior 
of the gluon condensate (G^) in an Abelian magnetic field is also nontrivial. Gluons do 
not carry an electric charge, but virtual quarks generated by them shift the quantity (G^) 
owing to their interaction with the magnetic field H. This phenomenon was studied in [35] 
making use of the low-energy theorems. The phase structure of QCD vacuum in an Abelian 
magnetic field at finite temperature was studied in [8, 15, 17, 19, 28, 29, 36, 37]. 

Relations arising from the symmetry properties play an important role in quantum field 
theories. Searches for symmetries and constraints that these symmetries impose on the 
physical properties of the system become particularly important in QCD as a theory that 
involves confinement, where composite states (hadrons) appear to be "observables". Low- 
energy theorems, or Ward identities (scale and chiral), play a fundamental role in obtaining 
deeper insight into the nonperturbative vacuum properties of QCD. In QCD, low-energy 
theorems were obtained in the early 1980s [38]. The low-energy theorems of QCD, which 
follow from general symmetry properties and which are independent of the details of the 
confinement mechanism, make it possible to obtain information that is sometimes inaccessible 
within any other method. Also, they can be used as "physically reasonable" constraints in 
constructing effective theories and various models of the QCD vacuum. In QCD, the low- 
energy theorems for T 7^ and Hq were obtained in [39, 40]. An important apphcation 
of low-energy theorems in a hot QCD was obtained in [41, 42]. Based on Kubo equation and 
low-energy theorems, the bulk viscosity of the quark-gluon matter was shown to be directly 
connected to a bilocal correlator of the energy-momentum tensor, and its value for a hot 
QCD case was evaluated [41, 42]. 

A method that makes it possible to generahze the low-energy theorems of QCD to the 
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case of finite temperature, a finite baryon chemical potential, and a nonzero magnetic field is 
developed in the present study. This method is used to study the nonperturbative vacuum 
and to derive an expression for a bifocal correlation function for the energy-momentum tensor 
and bulk viscosity in QCD for T 7^ 0, /x, 7^ and H 0. 



2. LOW-ENERGY THEOREMS OF QCD AT FINITE T, n AND H 

In the Euclidean formulation, the QCD partition function in presence of an external 
Abelian field can be represented in the form (here, T = 1/^ is temperature) 

Z = exp dxA J d^xF^^^ J [DB][Dq][Dq] exp |- ^ dx^ J d^xC^ , (1) 



where the QCD Lagrangian in a background field has the form 

q=u,d,... 

Here, Qg is the charge matrix for quarks of flavor q — {u, d,s,.. .) and bare mass moq and 
/iq is the quark chemical potential; the ghost and gauge-fixing terms are not written down 
explicitly in order to avoid encumbering the presentation. The pressure in the system (minus 
a thermodynamical potential) is determined by the expression /3VPo (T, fig, H, mog) = In Z. 
From the partition function in Eq. (1), one can obtain the following expression for the gluon 
condensate (G^) = {{Gl^f): 

dP 

{G'){T,^,g,H,mog) = -4g^y (3) 

The system described by the partition function in (1) is characterized by the set of di- 
mensional parameters Aq, //g, T, if, mog(Ao) and by the dimensionless charge ^IqC-^o), where 
the bare quark masses niog and the coupling constant are specified at the scale of the 
ultraviolet-cutoff mass Aq. On the other hand, we can go over to the renormalized (physical) 
pressure P and, with the aid of the properties of the renormalization group invariance of 
P, recast expression (3) into a form that involves derivatives with respect to the physical 
parameters T, /ig, and H and with respect to the renormalized masses nig. 

The dimensional transmutation phenomenon leads to the appearance of the nonpertur- 
bative dimensional parameter 

A = Aoexp|r -J^l, (4) 



as (A, 



0) 
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where Aq is the ultraviolet-cutofF mass, ag — qI/^t^ and ^{as) — dQ;s(Ao)/dln Aq is theGell- 
Mann-Low function. The quark mass mog has the anomalous dimension ^jyiq and depends 
on the scale Aq. The renormalization group equation for the running mass moy(Ao) has the 



form (ilnmog/rfln Aq = — Tm,, and we use the modified minimal-subtraction MS scheme, 
where j3 and are independent of the quark mass. The expression for the renormalization 
group invariant mass then has the form 

ruq ^ moq{K) Qy^V { j das > • (5) 

Since the physical (renormalized) pressure is a renormalization group invariant quantity, its 
anomalous dimension is zero. Thus, P has only a normal (canonical) dimension equal to 
four. Employing the renormalization group invariance of the quantity A, we can write P in 
the most general form as 

p = A^f(- N ^ (Q) 
•^^A' A' A2' A ^ ' 

where / is a function of the dimensionless ratios T/A, .... Prom (4), (5) and (6) we then 

obtain 

dP dP dA dP dmg 



d{l/gl) dA d{l/gi) ^ druq d{l/gl) ' 



Further, The anomaly in the trace of the energy-momentum tensor in QCD is related to the 
gluon condensate by the equation 

Taking into account (3), we obtain the gluon part of the trace of the energy-momentum 
tensor in the form 

K.) = [T§f + E".!; + 2// A + + _ 4) P. (10) 

Here and below, the energy-momentum tensor {9^^), the condensates (G^) and {qq), and the 
thermodynamic pressure P are functions of T, fig,H and m^. 

In the one-loop approximation, we have f3{as) — >■ —hal/2'K and 1 -|- 7^^^ — > 1, where 
h — (llA^c — 2A^/)/3. Thus, the gluon and quark parts of the trace of the energy-momentum 
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tensor in hot and dense QCD in a magnetic field can be expressed in terms of the physical 
pressure in the one-loop approximation as 

K,) -Y.^M) - -T.^gQ^ ■ (12) 

q q 1 

In the vacuum, that is, at T — 0, /ig — and H — 0, we arrive at the well-known 
expression for the nonperturbative vacuum energy density. In the one-loop approximation, 
this expression has the form 

Sva. = \{0l^ + Ol^)o = -P{T = 0, = 0, = 0, m,) 

--Y^2iG')o + \j2^q{qq)o. (13) 

Q 

By using the relations presented above, one can derive low-energy theorems of QCD at 
finite temperature and finite density in the presence of a magnetic field. Strictly speaking, 
the ^-function depends on if, so that the low-energy theorems could involve electromagnetic 
corrections, which are proportional to oc e^, but, since the physical pressure is independent of 
the scale Aq at which ultraviolet divergencies are regularized, one can choose an ultraviolet 
scale in such a way that Aq ^ H,T, fig, A. We can then restrict ourselves to the lowest 
order in the expansion of the /3-function, with the result that electromagnetic corrections 
disappear. Taking this into account, we consider the trace of the energy-momentum tensor 
in hot and dense QCD in the one-loop approximation, 

0,,{x) = -^(G^,(x))2 + Y,mqqq{x). (14) 
Also, we introduce the operator D, defining it as 

From relations (11) and (12), we find for the total vacuum expectation value of the trace of 
the energy-momentum tensor with allowance for massive quarks that 

{0,,) = {Ot, + 91^) = -^,{G') + Y.^M) = {D- 4)P. (16) 
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Differentiating Eq.(3) n times witfi respect to {I/Qq) and taking into account relations (6), 
(14), (15) and (16), we obtain 

= J dW..j ci^Xi(^^^(a:„)..J^^(xO^^,(0)). (17) 

To tfie rigfit-fiand side of (17) only connected diagrams are included, as usual. 

Similar arguments apply to an arbitrary operator 0{x) constructed from quark or gluon 
fields; that is, 

= j d'xn... j d'x,{ei^{Xn)...eUx,)dm, (18) 

where d is the canonical dimension of the operator O. If the operator O has an anomalous 
dimension as well, it is necessary to take into account the corresponding 7-function. 

Let us now consider the case of n = 1, which is of importance for physical applications. 
In other words, we will examine a bilocal correlation function for the tensors of the energy- 
momentum density in hot and dense QCD in a magnetic field. In terms of this correlation 
function, one can express the bulk viscosity of quark-gluon plasma in a magnetic field. 
Then, for gluon and quark contributions to the bilocal correlator of the trace of the energy- 
momentum tensor, we have the following relations: 

j d'x{ei{x)ei{o)) = {D- 4)(^^^), (19) 
j d'x{ei^{x)ei^m ^{D- 3){ei^). (20) 

Hence, for the the bilocal correlator of the trace of the energy-momentum tensor 

n = y" d^x{0^^{x)e^^io)) 

= / d'x{ei9f,^) + 2 I d^OM + 0{m% (21) 

where we included a correlator of quark summands to O(m^), being proportional to a quark 
mass squared; in what follows we will not take it into account. Based on relations (20), we 
find 

n = (^ - 4)(C) + 2(1) - 3)(0 -{D- A){e,,) + {b- 2){9V. (22) 
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3. BULK VISCOSITY ({T, /i, H) 

As it was shown in [42], according to a general Kubo formula, a bulk viscosity can be 
evaluated as a static limit of a bilocal correlator of the trace of the energy-momentum tensor. 

1 1 1'°" f 

^ ^ 9 i^o X ■ ^^^^ 

One can introduce a spectral density expressed in terms of a retarded Green function for the 
trace of the energy-momentum tensor 

p(a;,p) = --ImG^(c^,p). (24) 

TT 

Then, as was suggested in [42], for small frequencies the spectral density has the following 
form: 

P(^,0) 9C c^g .... 
= ~> — ^' (25) 

UJ TT CJq + <^ 

where the ouq parameter determines a scale at which a perturbation theory becomes appli- 
cable. Using this ansatz, the bulk viscosity can be written as 

oo 

9Ca;o - 2 J du ^ J d'x{e^,{x)e,^{0)) = U. (26) 



Thus, a problem to find the bulk viscosity ( reduces to a problem to evaluate the bilocal 
correlator 11. 

We extract from the correlator 11 a vacuum term. For this purpose, we write the following 
expression for the total pressure 

P^-e^ + P., (27) 

where is the pressure pure thermodynamical part. Quark and gluon contributions to the 
trace of the energy-momentum tensor can be written as 

W = (Oo + (28) 

and, using equation (16) 

(0,,) = + 91) = ie^ + {D- 4)n (29) 
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Allowing for the thermodynamic relation 

and taking into consideration the magnetic moment M — dP/dH we have 

{9^^) = 4£,,e + (£ - 3P), + 2MH. (31) 
Then the correlator (21, 22) can be written as 

n = no + n^ + n^, (32) 

where the vacuum contribution 

Ho = -4(^^^)o - 2{ei^)o = -16£vac - 2 5^m,(gg)o. (33) 
For the quark contribution HI we get the following expression 

The gluonic part Ilf of the correlator can be written as 



+ (r|^ + EN|; + 2i^J^-4)2Mi7. (35) 

We use the following definitions of thermodynamic quantities, in terms of pressure and 
energy density: entropy density s = dP/dT, specific heat = de/dT, velocity of sound 
= dP/de = s/cl, and medium magnetic susceptibility x = OM/dH = d^P/dH^. Then 
we find 

+ixH'-12MH + iH^T-^ + J2l',-^^M. (36) 

Thus, equations (33), (34), and (36) express the bilocal correlator of the trace of the energy- 
momentum tensor H, and, correspondingly, in accordance with (26), of the bulk viscosity ( 

via thermodynamic parameters: T, Hg, H, (gg)*, {e — 3P), s, c^, M, x 
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H{T, H)ujQ = -16£vac - 2 rnq{qq)o 

+ixH^-12MH + 4,H ^T-^ + Y,fi,-^jM. (37) 

We explore various limiting cases. Let us consider fig = H — and T 7^ 0, which 
corresponds to hot QCD studied in [42j. Taking into account (13) and using PCAC relation 

E ^Mq)o = -F^M'^ - FlMl (38) 
we write the vacuum contribution Ho as 

no = -4(e^^)o - 2(e^^)o = -4(e^^)o - 6(e^^)o = l^\eU + 6(F,2m,2 + M^), (39) 

where e^^^ is the gluonic contribution to the vacuum energy density (13). Assigning — 
if = in (37) and allowing for (38), (39), we find for the correlator 11 and , correspondingly, 
for the bulk viscosity 

9C(T)a;o ^Ts{^^-^- A{e - 3P% + (t-^ ^ 2^ Y^mq{qq). 

+16\elJ + &{FlMl + FlMi), (40) 

which exactly conforms to the main result in [42]. 

Now consider the utmost case of the cold quark matter T = H = and /Xg 7^ 0. And take 
into account that for such a variant 

^-3^= (>^/^<?^-4|i^ = y^/^,n,-4P, (41) 



where Ug = dP/dfiq is the quark density. We use relation (28) and definition P{^) = 
~£vac + P*, then we suppose T — H — in (37) , and, after simple transformations, we find^ 



^ The expression obtained for the bulk viscosity of the cold quark matter (42) does not coincide with a 
corresponding result in [43]. 
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9C{i^)uJo = 16P(/x) - 2 mq{qq) + ^ rriqi^q 



d{qq) 



dUgdUg, 



(42) 



We note that in (42) an expression for the all quark condensate {qq) including vacuum term 
is entering. 

In the limiting case of hot and dense quark-gluon matter, an expression for C(T, iJ,,H — 0) 
is obtained immediately from equation (37) for ({T,fj,,H) by substituting H = 0. 



According to Eq. (37), the bulk viscosity C iii ^ magnetic field is proportional to the 
magnetic susceptibility x- two different phases of strongly interacting matter (hadron 
and quark-gluon), the magnetic susceptibility has different signs [8]. This is because hadron 
matter at temperatures below the quark-hadron phase transition primarily consists of a gas 
of hot TT-mcsons, which in the magnetic field behaves as a diamagnetic medium of spinless 
charged particles. The system above the phase transition point is in the paramagnetic phase 
of hot quarks and gluons. This significantly affects the bulk viscosity of strongly interacting 
matter in the magnetic field, which should have a jump at the transition point, leading to 
interesting observable phenomena. 

In this paper a nonperturbative QCD vacuum at a finite temperature and a finite baryon 
chemical potential in an external magnetic field are considered. Relations between non- 
perturbative condensates and thermodynamical pressure at T 0, Hq ^ 0, and H ^ 
are obtained, and low-energy theorem are derived. A common formula for the bulk vis- 
cosity of the quark-gluon medium at T 7^ 0, /iq 7^ 0, and 7^ 0, is derived, which links 
C,{T,iJL,H) to thermodynamic system parameters: entropy density s — dP/dT, velocity of 
sound = dP/de — s/c^, non- ideality {e — 3P), quark condensate {qq)*, and medium mag- 
netic susceptibility x = dM/dH = d'^P/dH^. Some physically interesting limiting cases 
for the bulk viscosity are considered. These phenomena and exact relations obtained here 
may prove to be of importance in examining the quark-hadron phase transition in heavy ion 
collisions and in the early Universe, where strong magnetic fields are generated. 

I am grateful to A.B. Kaidalov, V.A. Novikov, and Yu.A. Simonov for useful discussions 
about the problems considered in the present study. 



4. CONCLUSIONS 
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